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Spin- Hall Effect and Diamagnetism of Dirac Electrons 
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Spin-Hall conductivity (SHC) of fully relativistic (4x4 matrix) Dirac electrons is studied based on the Kubo formula 
aiming at possible application to bismuth and bismuth-antimony alloys. It is found that there are two distinct contri- 
butions to SHC, one only from the states near the Fermi energy and the other from all the occupied states. The latter 
remains even in the insulating state, i.e., when the chemical potential lies in the band-gap, and turns to have the same 
dependences on the chemical potential as the orbital susceptibility (diamagnetism), a surprising fact. These results are 
applied to bismuth-antimony alloys and the doping dependence of the SHC is proposed. 
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The spin-orbit (SO) interaction has been known to lead 
to non-trivial and intriguing phenomena such as anisotropic 
magnetoresistance 1 ' 2 and anomalous Hall effect. 3 ' 4 More re- 
cent example is the spin-Hall effect (SHE) which indicates 
the generation of spin-current perpendicular to an external 
electric field. 5-11 One of the interesting features of SHE is 
the possibility of finite contributions even in the insulating 
states called as spin-Hall insulators. 9 Basically, materials with 
strong SO interaction should be the stage to observe such in- 
teresting features associated with SHE. In this context, bis- 
muth and its alloys with antimony will be one of the best can- 
didates since SO interactions are known to play crucial roles 
in these systems. 12 

Crystalline bismuth is basically cubic but with slight distor- 
tions (Peierls distortions) along the diagonal directions lead- 
ing to semi-metallic electronic states with small number of 
electrons (L-points) and holes (T-point). 13 Its electronic state 
has very characteristic features of anisotropy of small effec- 
tive masses and large g-factors (g ~ 100-1000), 14 ' 15 resulting 
from both small band-gap and strong SO interaction. Such 
characteristic electronic states are described by a 4 x 4 matrix 
Hamiltonian (conduction and valence bands with spin degrees 
of freedom) around L or T points in the k ■ p representation. 
This Hamiltonian, which may be called Wolff Hamiltonian, 
describes the low-energy properties of bismuth quite well, and 
turns out to be essentially the same as that for Dirac Hamilto- 
nian except with the anisotropy of velocity here, 16 because the 
space in solids is in general anisotropic in contrast to true vac- 
uum. (The Wolff Hamiltonian reduces to well-known Dirac 
Hamiltonian if spatial anisotropy of velocity is ignored, but 
with different velocity.) 

Bismuth is known for its very large diamagnetism with the 
particular feature of having maximum in the insulating states 
(by changing the chemical potential with alloying) whose ori- 
gin is now understood based on this Wolff Hamiltonian as due 
to inter-band effects of magnetic field in the presence of strong 
SO. 12 This inter-band effect also affects the off-diagonal trans- 
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port coefficient, but not the diagonal transport, in general. 17 
The inter-band effect in the off-diagonal transport coefficient 
and that in the diamagnetism should be related to each other 
in some way, but not clearly understood yet — a longstanding 
problem. 17 

In this paper, we study the SHE based on the Kubo formula 
for the isotropic Wolff Hamiltonian, 18 ' 19 i.e., the fully rela- 
tivistic Dirac Hamiltonian, not only in insulating but also in 
conducting states on an equal footing. We find that there are 
two distinct contributions to SHE, one from the states near 
the Fermi energy and the other from all the occupied state. 
The latter turns to have the same dependences on the chemi- 
cal potential as the orbital diamagnetism indicating that there 
is a close relationship between the SHE and diamagnetism in 
bismuth. 

We start from a simple one electron Hamiltonian with the 
SO interaction 

p 1 v 2 y l 

2m ^sm L c L 4m z c z 

as an effective model of electrons at L point in bismuth. 20 
Here the Pauli matrix cr corresponds to the real spin of elec- 
trons, V is the crystal potential, and the last term expresses the 
SO interaction. This effective Hamiltonian can be transformed 
into in an essentially identical form to the Dirac Hamiltonian 
as is shown by Wolff. 16 Here we discuss the isotropic case of 
the Wolff model: 1819 



\yk ■ cr 



-iyk ■ cr 



(2) 



where 2A is the band gap. The SO interaction is included in 
the matrix elements y, and the parabolic part in eq. (1) is dis- 
carded since it is negligibly small compared to the SO term 
in case of bismuth. 1415 Originally, the matrix elements are 
anisotropic, but we have assumed that all matrix elements are 
equal in order to make our arguments as simple and transpar- 
ent as possible in eq. (2). The eigen energy of this Hamilto- 
nian is +Ek = + ^y 2 k 2 + A 2 . The velocity operator is defined 
by v = dJt?/dk. 

As in the Dirac theory, the magnetic moment of electrons, 
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where g* = 2m e y 2 /A is the effective g-factor, m e is the free 
electron mass, and yUe = e/2m e c the Bohr magneton. Note that 
g* is given as the reciprocal of the effective cyclotron mass, 
g* = 2m e /m* c , where ra* = A/y 2 . It should be emphasized here 
that the sign of the magnetic moment are opposite between 
the conduction and valence band for the Dirac electrons. Cor- 
respondingly, the spin-velocity operator, v SI (i = x, y), which 
can be defined by the velocity of the magnetic moment along 
the z-direction, is given as 19 
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This is an Hermitian operator. 

The spin-Hall conductivity (SHC) is given as a linear re- 
sponse of the spin- velocity operator to the electric field on the 
basis of the Kubo formula: 



- [<D SVA (w + \S) - O SVA (0 + iS)] , 
CO L ' J 

ie„)v sy ^(ie„_)v x ] , 



r^Tr 

n.k 



(5) 
(6) 



where e„_ = s n - o>a, and ^(ie„) = (ie„ - is the Green 

function of eq. (2) with is„ = ie„ + irsgn(e„). Here we have 
introduced phenomenologically a quasiparticle damping rate 
T as an imaginary part of the self energy. 

After some straightforward calculations, we have 

4w e y 4 (s„ - £„_) 



<& m (\a) A ) = -eTY J -r 
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(7) 



We perform n-summation by the standard analytic continua- 
tion technique, and carry out the momentum integration. Then 
we obtain the following formula of the SHC 
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(9) 

(10) 



where e+ — e + iF, /(e) and /'(e) are the Fermi distribution 
function, whose energy is measured from the chemical poten- 
tial ju, and its derivative, respectively. (The branch cut of the 
square root is taken along the positive real axis.) Note that the 
diagonal spin conductivity, tr sxx is exactly zero in the case of 
the present Hamiltonian. 

In deriving (8)-(10), we have four contributions coming 
from the contours C\~\ in the complex z plane: C\ (Ci) is the 
contour from -oo to +oo (from +oo to -oo) along just above 
(below) the horizontal line Imz = cox, and C3 (C4) is the con- 
tour from -00 to +00 (from +00 to -00) along just above (be- 
low) the horizontal line Imz = 0. Kl vx in eq. (9) originates 
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Fig. 1. (Color online) Chemical potential dependence of the SHC (K^ )x ) 
for 17 A = 0.01 and E c / A = 100. The inset shows the plot of /f''" for different 
damping rates: T/A = 0.02,0.04,0.06,0.10. 



from the contribution of C2 + C, . It has a functional form 

1 r°° 

J d£f(£) [ GR(e)cA(e - w) _ ° R(e + ^ gA ^\ ' 

(id 

where G R(A) is the retarded (advanced) Green function. Con- 
fining ourselves to the static response, we need only w-linear 
term. Shift of the variable e — » e - co in the second term of 
eq. (11) leads to f JZ de {^f) G R (s)G A (s) in the limit of 
co — > 0. Therefore, the contribution of K\ yx is only from the 
states near the Fermi energy. This is similar to the contribu- 
tion in the transport properties, which may be called "trans- 
port contribution". 

The second term, Kf xx , in eq. (10) originates from the 
contribution of C\ + C4. It has a functional form of 
cof{s) [G R (s)G R (s) - G A (s)G A (s)]. This is similar to the con- 
tribution in the thermodynamical quantities such as thermody- 
namic potential, which may be called "thermodynamic contri- 
bution". In this case, the factor /(e) sums up the contributions 
from all the states below the Fermi energy. 

For the clean limit, F — > 0, at zero temperature, eqs. (9) and 
(10) can be expressed in the simple forms: 



' K lyx 
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(12) 



(13) 



where E c is the energy cutoff for the integration, and we dis- 
carded 0(A 2 / E 2 )-term. 

The chemical potential dependences of kL x and Kf vx are 
shown in Fig. 1 for which the values of T/A and E c /A are 
chosen to 0.01 and 100, respectively. K^ x gives the dominant 
contribution for the SHE. Furthermore, KR X , and so the total 
cr syx become maximum in the insulating region, \fi\ < A, and 
decays logarithmically for \fi\ > A although the carrier number 
increases. F-dependences of are shown in the inset of Fig. 
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1. The small F-dependence of both and K^ x (and so the 
total o~ syx ) means that the present SHE corresponds to the so- 
called intrinsic SHE. 711 

Now, let us discuss the inter-band effect on the SHE. 
We can divide Q^iiojj) into the contributions from the 
intra- or inter-band effect using the formula: 19 Q> S y X (io)A) = 
Z ajg <£lMi) with 



\k > A 



-eT^[(tfr a \v sy \(f/p)(^\v x \ilf a )%(ie n )%(is t ^)] , 

(14) 



n.k 



where a,B — + denote the conduction (+) or valence (-) band, 
t]/ ± are their wave functions, and ^±(is„) = [is„ + ■ This 
formula is equivalent to eq. (6). The inter-band contribution 
is from <Ifc and <&$y X , where the current and spin current ver- 
tices connect the conduction and valence bands. By explicit 
calculation of eq. (14), we find that the SHC originates only 
from the inter-band part, namely, 



o 
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(15) 
(16) 



where o$ x = + i5) - 0^.(0 + i5)] / (ico) . 

The unusual properties of <r m can be understood by com- 
paring with cr xx and cr vt 18 in terms of intra- and inter- 
band contributions. In the case of the diagonal conductivity, 
&xx(m> 0), the dominant contribution comes from the intra- 
band part, and it is suppressed as T increases (cr'"' ra oc F _1 ). In 
the case of the Hall conductivity, cr yx (a> — > 0), the dominant 
contribution also comes from the intra-band part, although 
there is a finite contribution from the inter-band part, which 
remains finite for T — > and exhibits anomalous properties 
similar to that of the orbital susceptibility. 18 The intra-band 
part is suppressed as cr 1 "" 3 oc r 2 , whereas the inter-band part 
is not. In the case of the SHC, which is completely inter-band 
effect, the F-dependence of <T syx (u> — > 0) is quite small as 
shown in Fig. 1 inset. From these results, we think that the 
T-dependence of the inter-band contribution will be generally 
small. 

Next, we discuss the relationship between cr syx and large 
diamagnetic susceptibility, x- Surprisingly, we find that the 
obtained has exactly the same /i-dependence x- In the 
limit of weak magnetic field, x is calculated by the simple but 
exact formula 21 

2 

x = ~ T X Tr [^(ie)v^(i£)v^(ie)v^(is)v y ] . (17) 
i,* 



Application of this formula to the present Hamiltonian results 
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(19) 



Despite their different starting points (eqs. (7) and (18)), the 
final expressions are equivalent; an astonishing result. As a 
consequence, we obtain the following relation between the 




x (Sb doping) 



Fig. 2. (Color online) Antimony doping, x, dependence of (a) cr syx (x) and 
— X i x ), an d (t>) electron bands at L-point and [i of bismuth. The form of A(,v) 
and n(x) are defined in the text. 



SHC and the orbital susceptibility in the insulating region : 

3f»eC 2 

OYvy = —y— X- (20) 

In the insulating region, cr xx is suppressed, so that there are 
no dissipative current. Only in such a case, the SHE becomes 
dissipationless, and it becomes exactly the same as the dissi- 
pationless diamagnetic current. We can say that the spin-Hall 
current under an electric field appears as the diamagnetic cur- 
rent under a magnetic field. This correspondence would be 
due to the duality between the electricity and magnetism. 

The relationship between the spin Hall conductivity and 
the spin density has been argued in the context of the two- 
dimensional quantum SHE of insulators for a non-relativistic 
model. 22 According to them, the "spin conserved" part of the 
spin Hall conductivity in the insulating case is given by a 
Stfeda-like formula: a\ff = -dS z /dB, where S z is the z- 
component spin density. Actually, in the present case, we find 
that the field derivative of the expectation value of magnetic 
moment /j ez (eq. (3)) is proportional to <x SVA as 



<r m oc -d{fis Z )/dB. 



(21) 



This is consistent with the claim of Ref. 22. However, (ji sz ) is 
not a measurable quantity and the relation in eq. (21) cannot 
be checked experimentally. Instead, eq. (20) can be examined 
experimentally. We also note that the Stfeda-like formula used 
in Ref. 22 is valid for the insulating case, while our calculation 
based on the Kubo formula is valid both for insulating and 
metallic case. 

Lastly, we discuss the implications of the present results to 
the experiments on bismuth. The band gap of electrons at L 
points is A = 7.7 meV, the chemical potential is ji — 35.3 
meV, so that ///A = 4.6 for pure semimetallic bismuth. 14 ' 15 
The band cut-off would be E c = 1-2 eV, 23 namely, E c /A = 
130-260. Substituting bismuth with antimony (Bi^SbJ can 
change the band structure as depicted in Fig. 2 (b). 24 By this 
substitution, the hole valence band at T point is lowered, re- 
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where the bonding (L s ) and antibonding (L a ) bands are in- 
verted, i.e., the topology of the bands changes from trivial to 
non-trivial. 25 ' 26 At x 0.07, the overlap between the elec- 
tron conduction bands and the hole valence band vanishes, so 
that <i reaches the band-edge, i.e., \i(x) = A(x). At x 0.09, 
the top of the hole valence band become lower than that of 
electron valence band, so that fi(x) = beyond this composi- 
tion. This doping dependence of /i and A can be simulated by 
±A(jc) = 1 - x/0.04 and = 4.6 - 4.6x/0.09. By substi- 
tuting these x-dependent A(x) and fi(x) into eqs. (8) and (19), 
the doping dependence of x(x) and cr syx (x) are obtained as is 
shown in Fig. 2 (a). 

The magnitude ofx(x) increases logarithmically, and a kink 
structure appears at x = 0.074 where fi reaches the band-edge. 
Beyond this composition, the magnitude of x(x) slightly de- 
creases due to the increase of A(x). These properties of x(x) 
agree well with the experimental results, 27 ' 28 indicating the 
validity of our theory. The magnitude of o- syx (x) exhibits a 
similar behavior to x(x); it increase logarithmically, and also 
have a kink at x = 0.074. There is no anomaly at around 
x = 0.04 at which the gap vanishes and the material is be- 
lieved to change from a simple insulator to a topological in- 
sulator. 25 ' 26 The difference between <x svc (x) and^f(x) is due to 
the intra-band contribution K* in cr syl . 

When we argue the transport property of bismuth, it is im- 
portant to consider the contributions from holes at T point, 
whose contributions are neglected in the present theory. The 
effective model for holes have been described by a parabolic 
band model with a large spin-mass term, 14 ' 15 which would be 
originated from the SO interaction. From this hole band, a fi- 
nite contribution is expected. However, this contribution will 
be small, since the gap at T point is much larger and the band 
cutoff would be much smaller than those of L points, as is the 
case of the hole contribution to the diamagnetism. 12 

In this paper, we have discussed the spin-Hall effect of the 
isotropic Wolff Hamiltonian, which is fully relativistic, on the 
basis of the Kubo formula. It has been shown that the spin- 
Hall effect appears in this Dirac electron system. Especially, 
the spin-Hall effect becomes maximum in the insulating re- 
gion, where the electric current hardly flows. Therefore, we 
can obtain the dissipationless spin-Hall current there. This 
spin-Hall conductivity originates only from the inter-band ef- 
fect, while the Hall conductivity originate both from intra- 
(dominant) and inter-band (small) effects. 

It has been found that there are two distinct contributions, 
i.e., K\ , "transport contribution", and Kf xx , "thermodynamic 
contribution", to the SHE. The latter has exactly the same 
chemical-potential dependence as the orbital susceptibility 
(diamagnetism). Eqs. (10), (19) and (20) are the first proof 
that shows the definite relation between the spin Hall conduc- 
tivity and the diamagnetism with the fully relativistic model. 
This astonishing correspondence strongly suggests the spin- 
Hall effect has the same nature as the diamagnetism through 
the inter-band effect associated with the spin-orbit interac- 
tion and the duality between electricity and magnetism. This 
would play a crucial role for understanding the long-standing 



problem: how the transport coefficients relate to the orbital 
susceptibility. 17 

The present work gives a fully relativistic theory of the 
spin-Hall effect not only in insulating but also in conducting 
state on equal footing for the first time. Also, our results have 

disclosed that the Dirac electron system will provide another 
ideal situation for investigating the spin-Hall effect, since it is 
very simple and gives not only qualitative but also quantitative 
agreements with experiments on bismuth. 
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